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Abstract 

A condensed introduction to quantum gauge theories is given in the perturbative 
S-matrix framework; path integral methods are used nowhere. This approach empha- 
sizes the fact that it is not necessary to start from classical gauge theories which are 
then subject to quantization, but it is also possible to recover the classical group struc- 
ture and coupling properties from purely quantum mechanical principles. As a main 
tool we use a free field version of the Becchi-Rouet-Stora-Tyutin gauge transformation, 
which contains no interaction terms related to a coupling constant. This free gauge 
transformation can be formulated in an analogous way for quantum electrodynamics, 
Yang-Mills theories with massless or massive gauge bosons and quantum gravity. 
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1 Introduction 

It is well known that gauge theories play a fundamental role in modern quantum field 
theory. Prom this observation one might conclude that gauge invariance is a physical 
mechanism inherent in many particle interactions. But this is not the case. Gauge invari- 
ance is rather a mathematical artefact which stems from the way we formulate (quantum) 
field theory, although its presence is very helpful in order to achieve a consistent formula- 
tion of many quantum field theories. Like in the case of general relativity, where we are free 
to choose the coordinate system according to our taste, we also have some freedom how 
to deal with field theories. The free quantum fields which are used in perturbation theory 
in order to calculate physical quantities like e.g. cross sections simply provide a 'coordi- 
natization' of the problem under consideration, and they are 'as physical' as a coordinate 
system. 

One may illustrate this by a simple observation. Massless particles with spin have 
only two degrees of freedom, e.g. photons appear only in two helicity states. Such states 
are called 'physical' in the paper, whereas timelike or longitudinal photons, which are 
an unavoidable byproduct of covariant quantization, are called 'unphysical'. But if we 
take into account that one can introduce running coupling constants and running masses 
in interacting field theories, it becomes clear that also 'physical' particles like quarks are 
merely mathematical constructions, since we are always free to change the structure of the 
Fock space underlying our calculations by renormalizing the mass of the quark asymptotic 
states. Therefore it is meaningless to ask how a composite particle can be decomposed 
into free 'naked' particle states with a well-defined mass. What one really should do is 
to calculate vacuum expectation values of products of local interacting fields in a nonper- 
turbative manner. From these distributions it is in principle possible to reconstruct the 
true physical Hilbert space, as it is claimed by the reconstruction theorem of Wightman 
d]. But this is a very difficult task and has been performed only for exotic cases like e.g. 
quantum field theories in 1-1-1 spacetime dimensions [H 01 13 EI- 

It is not the intention of this introduction to give a complete description of all aspects 
related to the arbitraryness of the formulation of quantum field theory, and for mathe- 
matical details we must refer to the literature. But since we are working in a strictly 
perturbative sense by using free field operators only, the whole formalism used in this 
paper can be formulated in a fully mathematical manner if necessary. 

The structure of this paper is as follows: First, we start with the free classical elec- 
tromagnetic field (or, more precisely, with the classical gauge theory of a massless spin 
1 field), for which the basic issues of classical gauge invariance are explained. Then the 
electromagnetic field is subject to quantization, and we lift the gauge transformation on 
the quantum level. It is explained how the quantum gauge transformation, which is the 
free field analogue of the full Becchi-Rouet-Stora-Tyutin (BRST) transformation [HI EI; 
can be applied to quantum chromodynamics (QCD) along the same lines as for the quan- 
tum electrodynamics (QED) case. As a further step, the method is extended to massive 
gauge theories, where the Higgs field is involved. Finally, we outline how quantum gauge 
invariance can be formulated for quantum gravity. 
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2 The classical electromagnetic field 

The equation of motion for the noninteracting vector potential A^ in electrodynamics can 
be derived from the Lagrangian 

C = -^^F,,F^'' = ^{E'-B% (1) 

where F^u = d^A^ — d^A^ is the electromagnetic field strength tensor which can be 
expressed directly in terms of the components of the electric and magnetic fields E and 
B. From the Euler-Lagrange equations 

one derives then the wave equation 

d.d-'A^ - d^d.A^ = 0. (3) 

A problem arises from the fact that the physically measurable electromagnetic fields E and 
B remain unchanged if the vector potential is subject to a gauge transformation expressed 
by a real scalar field u {x = {x^ , x)) 

A^{x) ^ A^{x) + di,u{x), (4) 

i.e. there is a redundance in the local description of the electromagnetic fields by a vector 
potential, aside from global topological information which may be contained in the poten- 
tial like in the case described by Aharonov and Bohm 0. From the mathematical point 
of view, we generally assume that the objects we are dealing with 'behave well', i.e. A^ 
and u should be differentiable and behave at large spatial distances in such a way that 
physical quantities like the total energy of the field configuration are finite. 

The possibility to transform the vector potential without changing its physical content 
allows to impose gauge conditions. Due to its manifestly covariant form, a very popular 
choice is the Lorentz gauge condition, d^A^ = 0. This condition can be enforced by 
transforming the initial vector potential according to @ with a field u which fulfills 

□n = d^.d^'u = -d^A^". (5) 

Then the new transformed field satisfies the simple wave equation 

□A'' = d^d'^A'' = 0. (6) 

But the vector potential is not completely determined by the Lorentz gauge condition. 
It is still possible to apply a gauge transformation to the vector potential if u fulfills the 
wave equation d^d'^u = 0, such that the transformed vector potential is still in Lorentz 
gauge. In the sequel we will denote the field n, which is somehow related to the unphysical 
degrees of freedom of the vector potential in Lorentz gauge, as 'ghost field'. The reason 
will become clear at a later stage. 
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It is important to note that the equation of motion for the vector potential can be 
manipulated by a modification of the Lagrangian according to an early observation by 
Fermi |^ . Adding a gauge fixing term to the Lagrangian 

C = -\f,,F^'^ - lc{d,An', (7) 

where is an arbitrary real parameter, leads to the equation of motion 

nA^ - (1 - Odf.d.A" = 0. (8) 

One may argue that the gauge fixing does not change the physical content of the La- 
grangian, since we may enforce the Lorentz gauge condition df^A^ = 0, and then the 
gauge fixing term vanishes in the Lagrangian. For the special choice ^ = 1 we recover 
the wave equation ©. This choice is often referred to as Feynman gauge, although this 
terminology is a bit misleading. We must clearly distinguish between the gauge fixing, 
which fixes the Lagrangian and determines the form of the wave equation (e.g. the Feyn- 
man gauge), and the choice of a gauge condition for the vector potential (e.g. the Lorentz 
gauge), which reduces the redundance originating from the use of a vector potential as 
a basic field. It is possible to impose more general (nonlinear) gauge fixing conditions, 
e.g. by adding a term ~ {d^A^^Y to the Lagrangian, but we restrict ourselves to the most 
important cases here. We conclude this discussion of the gauge fixing with the remark 
that it is indeed necessary to add appropriate gauge fixing terms to the free Lagrangian 
(PQ), since otherwise problems arise in the construction of the free photon propagator. 

Up to now, the discussion of the electromagnetic field has been completely classical. In 
the next section, we will consider a quantized version of the free electromagnetic field. For 
the field operator A^{x) we will use the Feynman gauge such that the operator obeys the 
simplest equation of motion ©. All calculations could be performed in a strictly analogous 
way for arbitrary ^—gauges, but the choice ^ = 1 has obviously notational advantages [TU] . 
Strictly speaking, each choice of ^ leads after quantization to a different quantum field 
theory with different Fock spaces and field operators. But physical observables like cross 
sections etc. should be independent of the gauge fixing. 

3 Quantization of free fields 

It is well known that a free massless hermitian scalar field ^{x) can be represented by 

ip{t, x) = (27r)-3/2 I (a(fc)e-'("*-^"^) + a(fc)^e'("*-^"^)) , (9) 



/2uj 

where w = and the creation and annihilation operators fulfill the commutation relations 

[a{k),a{k'y\ = 5'^'^\k-k'), (10) 

[a{k), a{k')] = [a(A?)t, a(P)t] = 0. (11) 

The ^ denotes the adjoint with respect to a positive definite scalar product so that the 
operators can be represented in the usual way in a Fock space with unique vacuum, 
is the three dimensional Dirac distribution. 
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We now try to quantize the vector potential A^(x) as four independent real scalar 
fields 

Af'it, x) = (27r)-3/2 f ^ (^a^(A?)e-^(^*-^^) + a'^(fc)te*('^*-^^)) , (12) 



'2u; 

and naively we assume that the commutation relations for the creation and annihilation 
operators are 

[a^(A?), a^ik')] = [a^{k)la''{k')^] = 0. (13) 



= is automatically satisfied according to (|12|) . since the spacetime dependence 
of the field operator is given by the plane wave terms e^**^^ = e^*('^*~'^^), and 6 is the 
Kronecker delta of the indices fi, v. 

A calculation shows that the free massless hermitian scalar field ip fulfills the commu- 
tation relation 

[ip{x),^{y)]=-iD{x-y), (14) 
where D(x — y) is the (massless) Pauli-Jordan distribution 

D{x) = I d4fc5(A;2)sgn(fcO)e-*^^ = ^sgn{x^)6{x^). (15) 

The Pauli-Jordan distribution fulfills the distributional identity 

doD{x)l,,=o = 6'^^\x) (16) 

which implies the equal time commutation relation for a scalar field (p and its canonical 
momentum tt = (p 

Mx),^{y)]\„o = iS^''Hx-y) (IT) 

which is often taken as the starting point for the quantization of the scalar field. The 
distribution D{x — y) was introduced by Pauli and Jordan in 1928, when they performed 
the first covariant quantization of the radiation field . The naive choice (|12|) comprises 
a problem, since the commutation relations 

[Af'ix),A''iy)] = -i6^,D{x-y) (18) 

are not covariant. A simple way to remedy this defect is to change the definition of A^ 
into 

A^ix) = (2^)-3/2 f (aO(fc)e-^^" - a\k)^e''''') , (19) 



i.e. we make A^ a skew-adjoint operator instead of self-adjoint. Another strategy is to 
introduce a Fock space with negative norm states, as proposed by Gupta and Bleuler 
|121 113j . The commutation relations then become 

[A^{x),A''{y)]=ig^^''D{x-y). (20) 

We have introduced a nonhermitian field in order to save the Lorentz symmetry of the 
theory, and therefore one might expect problems with the unitarity of the QED S-matrix. 
But no problems arise, a fact which is related to the gauge symmetry of the theory. 
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A concluding remark is in order here. It is possible to avoid unphysical polarization 
states by quantizing the photon field in radiation gauge, which is expressed classically by 
the conditions 

A^ = 0, V-A = 0. (21) 

The (quantized) photon field operator contains then only physical (transverse) polariza- 
tions 




with 

e{k, A) • e{k, A') = 6xx' and e{k, X) ■ k = (23) 

and standard commutation relations for a and . But this is only seemingly an advantage 
because manifest Lorentz symmetry is lost in H22|) . and no general strategy is known how to 
deal with the renormalization of divergent higher order contributions in interacting theories 
without the guiding help of Lorentz symmetry. Thus we have to make the choice between 
having unphysical particles, but also manifest Lorentz covariance in our calculations, or 
working only on a physical Fock space, but without having a consistent scheme at hand 
to regularize loop diagrams in perturbative calculations. 



4 A simple version of a quantum gauge transformation 

We introduce now a quantized version of the classical gauge transformation for the free 
vector potential 

^ yl^ + df'u, nu = 0. (24) 

Of course there exist different approaches for the treatment of gauge invariance in quantum 
field theory, but the approach presented here is quite simple and does the job. 
As a first step, we define the gauge transformation operator or gauge charge Q 

Q= j cPxd^A^'{x)dQu{x). (25) 

XQ=const. 

It is sufficient for the moment to consider n as a real C-number field. In the case of QCD, 
u will necessarily become a fermionic scalar field. It can be shown that Q is a well-defined 
operator on the Fock space generated by the creation and annihilation operators according 
to (|13() . It is not important over which spacelike plane the integral in ()25() is taken, since 
Q is time independent: 

Q= j {-dld^A^'u + ^^A^'^lu) 

xo=const. 

(fx (- A ^^A^'u + d^A^" Au) = 0. (26) 

xo=const. 
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This formal proof uses the wave equation and partial integration. Another way to show 
the time independence of the gauge charge is to define the gauge current 

fg=d,A^d^u, Q = ld'xfg, (27) 

which is conserved 

= d^id^A'^d'^u - d'^d.A^u) = 0. (28) 

A crucial property of the gauge charge is expressed by the commutator with A^^, which is 
a C-number 

[Q,A^'{x)]=i^^'u{x), (29) 

and all higher commutators like 

[Q,[Q,^^(x)]] =0. (30) 

vanish for obvious reasons. The reader is invited to check the commutation relation H29|l . 
an outline of the calculation can be found in Appendix A. Therefore we have 

g-iAQ^Mg+iAQ = _ ^[Q,^^] - [Q,A^']] + ... 

= A^'- i\[Q, A^'] = A^' + Xd^'u, (31) 
i.e. Q is a generator of gauge transformations. 

5 Definition of perturbative quantum gauge invariance 

We take the next step towards full QED and couple photons to electrons. In perturbative 
QED, the S-matrix is expanded as a power series in the coupling constant e. At first order, 
the interaction is described by the normally ordered product of free fields 

Wint(^) = : ^{x)r^{x) : A^{x), (32) 

where ^ is the electron field operator. The S-matrix is then usually given in the literature 
by the formal expression (T denotes time ordering) 

S = i + Y.^- / d^Xi...d^XnT[n;^^{xi)...n;^^{Xn)] 
n=l ^' 



If 

1 + ^ ^ / d^Xl . . . d^Xn Tr 
n=l ' 



n(xi, . . . Xn), (33) 



where we have introduced the time-ordered products r„ for notational simplicity, and we 
have 

T,{x) = -m^nti^) = ie : ^{xh'^-i'ix) : A^{x). (34) 

Expression (|33|) is plagued by infrared and ultraviolet divergences (see Appendix E). We 
leave this technical problem aside and we assume that the r„ are regularized, well-defined 
operator valued distributions, which are symmetric in the space coordinates {xi, . . . Xn)- 
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In the previous section, we used a real C-number field u in order to show how the gauge 
transformation can be lifted on the operator level. We assume now that n is a fermionic 
scalar field, and additionally we introduce an anti-ghost field n, which is not the hermitian 
conjugate of u. The reason why we choose u, u to be fermionic, against all conventional 
wisdom that particles without spin are bosons, is very simple: Bosonic ghosts do not allow 
to formulate a consistent theory in the case of quantum chromodynamics, and therefore 
we start directly with the 'correct' strategy. A second reason is that we want Q to be 
nilpotent; = which allows the definition of the physical Hilbert space as 

^phys = ker Ql ran Q (35) 

(see Appendix C). 

We give now a precise definition of perturbative quantum gauge invariance for QED, 
which will work in a completely analogous way for QCD. 
Let 

Q -= j d^xd^.A^'ix) do u{x) (36) 

be the generator of (free field) gauge transformations, called gauge charge for brevity. This 
Q has first been introduced in a famous paper by Kugo and Ojima jl4|. The positive and 
negative frequency parts of the free fields satisfy the {anti-} commutation relations 

[4±)(x), AW(y)] = ig,,D^^Hx " v), (37) 



{nW(x),nW(y)} = -iD^^\x-y), (38) 

and all other commutators vanish. Here, are the positive and negative frequency parts 
of the massless Pauli-Jordan distribution with Fourier transforms 

Z)(±)(p) = ±^e(V)5(p'). (39) 

All the fields fulfill the Klein-Gordon equation with zero mass, and as already mentioned, 
we are working in Feynman gauge, but the following discussion would go through with 
some technical changes in other covariant ^-gauges as well. Note that the quantum field 
A^ does not satisfy d^A^^ = on the whole Fock space, but only on the physical subspace. 



An explicit representation of the ghost fields is given by 

u{x) = (2^)-3/2 J ^(^C2{k)e-'''' + cl(fe)e*^"), (40) 

u{x) = (2vr)-3/2| ^(^-c,{k)e-^''^ + cUk)e^'^), (41) 

{c^{k),c]{k')} = 6i,6'^^\k-k'), i,j = 1,2, (42) 



In order to see how the infinitesimal gauge transformation acts on the free fields, we 
have to calculate the commutators (see |^ and Appendix A) 

[Q,A^]=idf,u, {Q,u} = 0, {Q,u} = -id^A'', [Q, ^r] = [Q, ^] = 0. (43) 
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The commutators of Q with the electron field are of course trivial, since the operators act 
on different Fock space sectors. We need only the first and the last commutator in ()43() 
here, the others will become important in the QCD case. From we know that the 
commutator of Q with an operator gives the first order variation of the operator subject 
to a gauge transformation. Then we have for the first order interaction Ti 

[Q,Ti{x)] = -e : ^'(x)7^^(x) : a^n(x) 

= idf,{ie : ^{x)-r^''^{x) : u{x)) = idf,T^^^{x). (44) 
Here we used the fact that the electron current is conserved 

df, : ^'7'^^' := 0. (45) 

because ^ fulfills the free Dirac equation. Note that the free electron field is not affected 
by the gauge transformation. We may can r^*/^ = ie : ^'7^^' : u the 'Q-vertex' of QED. 
The generalization of ()44|) to n-th order is 

n 

[Q,Tn{xi, ...Xn)] = i'^d^'T^^iixi, ...Xn) = (sum of divergences) , (46) 
1=1 

where T^^^ is again a mathematically well-defined version of the time-ordered product 

r^/^(xi,...,x„)" ="r(ri(xi)...rf/,(x0...ri(x„)). (47) 

It is relatively easy to see how (|46p comes about if one understands the second order case. 
This example is discussed in detail in Appendix D. We define to be the condition of 
gauge invariance [TCj . 

If we consider for a fixed xi all terms in T„, with the external field operator Afj^{xi) 

Tn{xi, ...Xn) =■■ tf (xi, ...Xn)A^{xi) : +... (48) 

(the dots represent terms without Afj_{xi)), then gauge invariance (|I6|) requires 

^^[if (Xi, ...XnHxi)] = tf (Xi, ...Xn)d^u{xi) (49) 

or 

dltf{xu...xn) = 0, (50) 

i.e. we obtain the Ward-Takahashi identities ^H] for QED. The Ward-Takahashi identities 
express the implications of gauge invariance of QED, which is defined here on the operator 
level, by C-number identities for Green's distributions. 

We have found the following important property of QED: There exists a symmetry 
transformation generated by the gauge charge Q, which leaves the S-matrix elements 
invariant, since the gauge transformation only adds divergences in the analytic sense to 
the S-matrix expansion which vanish after integration ov6r th.6 coordiricitGS xi, . . . (see 
Appendix E). 
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The S-matrix 'lives' on a Fock space J- which contains physical and unphysical states. 
E.g., unphysical single photon states can be created by acting on the perturbative vacuum 
\Q) with timelike or longitudinal creation operators: 

a\k)^\n), kja^k)^\n), (51) 

i=l,2,3 

and states which contain ghosts are always unphysical. The physical subspace J^phys is the 
space of states \^) without timelike and longitudinal photons and ghost, that means 

a^(k)\^) = 0, kja^{k)\^) =0, ci,2(A?)|«>) = (52) 

Note that the definition of the physical space is not Lorentz invariant, but the definition 
(i35|) is. 

The observation that QED is gauge invariant is interesting, but the true importance 
of gauge invariance is the fact that it allows to prove on a formal level the unitarity of 
the S-matrix on the physical subspace (see the last paper of ^S])- Due to the presence 
of the skew-adjoint operator or the presence of unphysical longitudinal and timelike 
photon states, the S-matrix is not unitary on the full Fock space, but it is on Tphys- We 
do not give the algebraic proof here, but we point out that gauge invariance is the basic 
prerequisite which ensures unitarity. A detailed discussion of this fact can be found in 
|151ll7l[TK| . We have introduced the ghosts only as a formal tool, since they 'blow up' the 
Fock space unnecessarily, and they do not interact with the electrons and photons. But 
in QCD, the situation is not so trivial. 



6 Quantum chromodynamics 

As a further step, we consider QCD without fermions (i.e. quarks) now. The first order 
coupling of the gluons that we obtain from the classical QCD Lagrangian is given by (see 
Appendix B) 

T^ix) = i\fabc ■ A^a{x)A,,{x)F:^^'{x) : 

= igfabc:A^a{x)A^b{x)d''A^{x):, a,b,c = l,...8 (53) 

where Fc'^ = d^A^ — Ac and fabc are the totally antisymmetric structure constants of 
the gauge group SU(3). The asymptotic free fields satisfy the commutation relations 

[Afa^{x),A^X'{y)]=i5a,g,,D^^\x-y) (54) 

and 

{u^^\x)M^\y)] = -i5atD<^^\x-y), (55) 
and all other {anti-}commutators vanish. Defining the gauge charge as in (|36|) by 

Q:= j d^xd^A^ix) do ua{x), (56) 

where summation over repeated indices is understood, we are led to the following commu- 
tators with the fields: 



[Q,A^]=idf'ua, [Q,Fr]=0, {Q,M = 0, {Q,ua} = -idf^A^a ■ (57) 
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For the commutator of Q with we obtain 

[Q, T,^{x)] = igfabc : {[Q, A^^jA.bd-'A^ + A^a[Q, A.^WA'^ + A^^A.biQ, d^A^]} : 

= -gfabc : {d^UaA^bd'^A'^ + d^iA^aUbd^A^)} : . (58) 

The last term is a divergence, but the first term spoils gauge invariance and therefore the 
unitarity of the theory. To restore gauge invariance, we must somehow compensate the 
first term in H58() . 

We consider therefore the gluon-ghost coupling term 

T^{x) = igfabc : A^aix)ubix)df'ucix) :, (59) 

For the commutator with Q we get 

[Q, T^] = ig : [Q, A^a]ubd^'u, + A^a{Q, Ub}d^'u, - A^aUbd^'iQ, nj : 

= -gfabc ■■ {d^.UaUbd^'Uc + A^aUbOf^d" A^^} : • (60) 

Note that we can always decompose the commutator in a clever way such that only an- 
ticommutators of Q with the ghost fields and commutators with the gauge fields appear. 
Taking the antisymmetry of fabc and : UaUb := — : UbUa ■ into account, we see that the 
first term is a divergence 

fabc ■ dfj.UaUbd'^Uc := ^fabcdf, : UaUbd'^Uc : (61) 

because udx) fulfills the wave equation. The second term can be written as 

-gfabc ■■ A^,aUbd^d^ A,, := -gfabc ■ {A^^^Ubd^ A,,) - d'^UbA^^d^A,, : . (62) 
Interchanging a h and in the last term, it becomes equal to the first term in 

Hence, the combination 

Ti{x) = igfabc{\ : A^a{x)A,b{x)F^^'{x) : - : A^a{x)ub{x)d^'u,{x) :} (63) 
leads to a gauge invariant first order coupling 

[Q, Ti] = gfabc : -d^iA^aUbid^A^ - 9'^^^) + ^d^iuaUbd^Uc) : . (64) 

In contrast to QED, the nonlinear selfinteraction of the unphysical degrees of freedom 
of the gluons spoils gauge invariance. Coupling ghosts in an appropriate way to the 
gluons restores the consistency of the theory. Thats why we need ghosts. They interfere 
destructively with the unphysical gluons and save the unitarity of the theory. 

The first order coupling given in H63|) is in fact not the most general one, and it is 
also possible to construct a gauge invariant first order coupling with bosonic ghosts |15j : 
but gauge invariance then breaks down at second order of perturbation theory. Gauge 
invariance at order n is again defined by (|46j) . We have discussed gauge invariance only at 
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first order here. To prove that gauge invariance (and other symmetries of the theory) is 
not broken by renormalization, is the 'hard problem' of renormalization theory jl9ll2L)l[^ . 

The 4-gluon term ~ g'^ which also appears in the classical Lagrangian is missing in 
Ti. This term appears as a necessary local normalization term at second order, and 
its structure is also fixed by perturbative gauge invariance. This stresses the fact that 
perturbative gauge invariance is strongly related to the formal expansion of the theory in 
powers of the coupling constant g. 

We started this section by presupposing that the coupling of the gluons is already given 
by the classical Lagrangian. But perturbative gauge invariance is indeed a very restrictive 
condition: It fixes the interaction to a large extent. We outline here how this fact can be 
derived. Only a part of the derivation is given here, and it is not assumed that the reader 
will check the calculations in detail. A full discussion is given in UH], and for further 
reading we recommend also |22j . 

We start by describing the interaction of the massless gluons by the most general 
renormalizable ansatz (i.e. the dimension of the interaction terms is energy^) with zero 
ghost number (coupling with non-zero ghost number would affect the theory only on the 
unphysical sector) 

fiix) =ig{ f^,,:A^aix)A,,ix)d''A>^{x): + 

fibc ■■ Af^aUbd'^Uc ■■ + 
fabc ■■ A^adf'UbUc ■ + 
fabc '■ ■^fiaA'^duA^^ : + 

fabc '■ duA'^UbUc :| , fabc = fbac^ 

where the /'s are arbitrary real constants. This first order coupling term is then antisym- 
metric with respect to the conjugation K defined in Appendix C. 

Adding divergence terms to Ti will not change the physics of the theory, since diver- 
gences get integrated out in ^\VA\i . and we assume that this property is not destroyed by the 
renormalization of the theory. Furthermore, the (anti-)commutator of Q with a divergence 
is also a divergence. We can therefore always calculate modulo divergences in the follow- 
ing. We modify f, by adding -f (/i^, + f^Jd, : A^^A.bA^ : and -igf^bcd,. ■ A^^UbUc : 
and arrive at a more compact, equivalent first order coupling 

Ti =ig{ flbc-\aAubd'Ar.+ 

fabc '■ A^aUbd^Uc : + 
fabc '■ ■AfJ.aA^dyA'^^ : -|- 

flbc-duKubUc-], (65) 

where 

fabc = ~fcba fabc = fbac' (66) 
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The gauge variation of Ti is 

[Q,Ti] =-9fl,, {■.d^UaA,,d-A>^,+ 



9 fabc 



-^9ftbc 

^9fabc 



A^aduUbd'A'^ + 
Af^^A^bd-'d'^Uc :} 

I : d^UaUbd^Uc + 

A^aUbd^d.A^, :} 
: d^UaA^^d^A';^ : + 



+ 



+ 



(67) 



Gauge invariance requires that the gauge variation be a divergence, hence we make again 
a general ansatz 

[Q,ri] =gd^{ gl,, : d'^UaA.bA'', : + 

9abc • '^aA'^dyA'^^ : + 
9lbc ■ d^UaUbUc : + 

5^bc '■ UaUbd^Uc ■■ + 
9abc '■ duUaAlA'^ : + 
9abc '• Uad^ A\^Ayc '■ + 

where g\^^^ = g\(.b^9abc ~ ~9bac- Comparing the terms in (|67j) and (|68|) . we obtain 
constraints for the couphng coefficients: 



(68) 
a set of 



- flab = '^9lbc+9lbc from 


d'^udf.A^A'' 


(69) 


2,5 9 fi 


d^uAf.duA" 


(70) 


„2 , 2 _ .5 

ilabc ' ilacb J bac J cab 


ud^A'^d^A^' 


(71) 


~fabc = 9abc + ^^afec 


d^uud^u : 


(72) 


n3 - 
iJabc ilbac 


d^ud^uu : 


(73) 


_ f 1 _ f 1 - „5 , 5 
J cfea J bca ilabc ' Wacfe 


d^d^uA'^A^ 


(74) 


f 1 _ „5 , 7 
./ acb iJabc ~r i/afec 


df.ud^Af'A'' 


(75) 




ud^A^d^'A" 


(76) 


/ cab 9acb ' 9abc 


ud^.duA^'A'' 


(77) 


7 _ 7 
5abc fl'acfe 


ud^A^d^A"" 


(78) 



Prom ((711), IZSl), and we readily derive 

f abc ~^ f acb f bca f cba 



(79) 



Combining this with (|66j) . we obtain the first important result that /^j,^ is totally anti- 
symmetric. g\f^^ is symmetric in h and c; from (|76() and ()69() we conclude (/^^^ = and 



9abc 
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We did not yet fully exploit gauge invariance. Taking into account additionally that 
the gauge charge is nilpotent = and consequently 

{Q,[Q,ri]} = (80) 

provides further relations which fix the interaction to 

ri = rf^Vrf, (8i) 



Tl^ = igfLb ^A^aA^tFr : - : A^aUbd'^u, :}, (82) 

tF = igf!bc ■■ duAytUc :, (83) 
where 'YM' stands for 'Yang- Mills' and 'D' for 'Deformation', and /^^^ = — 

We stop our analysis here for the sake of brevity. A detailed analysis of gauge invariance 
at second order shows |18| that the coupling coefficients fabc uiust fulfill the Jacobi identity 

fabcfdec ~l" fadcfebc ~l~ faecfbdc ~ 0; (^^) 

i.e. we are lead to the nice result that the fabc are structure constants of a Lie group. 
Note that the Lie structure of gauge theories is usually presupposed, but here it follows 
from basic symmetries of quantum field theory, namely quantum gauge invariance (which 
is related to unitarity) and Lorentz covariance. 
Finally, we are left with the term 

Tf = igfl,^ : 9.A>6ne : . (85) 

One can show that this interaction term which contains only unphysical fields does not 
contribute to the S-matrix on the physical sector. 

The reason why gauge invariance at higher orders of the perturbation expansion is 
a nontrivial task can be understood qualitatively on a basic level. Commutators of free 
fields can be expressed by (derivatives of) Pauli-Jordan distributions, e.g. we have 

{d,A^{x),d^A^^{y)] = ig^^dldlD{x - y) 



= -in,^D{x - y) = 0. (86) 

But through the timeordering process, Feynman type propagators ~ Dp{x — y) appear 
in the amplitudes instead of Pauli-Jordan distributions, which fulfill the inhomogeneous 
wave equation 

□Df(z) = -5(^)(z), (87) 

i.e. derivatives acting on Feynman propagators may generate 'anomalous' local terms, 
because time ordering does not commute with analytic derivation. Therefore, maintaining 
gauge invariance may imply restrictions to the coupling structure also at higher orders. If 
it is not possible to fix the theory such that gauge invariance can be maintained, the theory 
is called nonrenormalizable. The result obtained by the analysis above can be considered 
as an inversion of 't Hooft's famous result that Yang-Mills theories are renormalizable |23j . 

In this section, we focused on the purely gluonic part of QCD. A full discussion would 
include also the coupling of gluons to quarks, but all important features are contained in 
the present discussion. 
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7 Yang-Mills theories with massive gauge fields 

In this section we give a short discussion of a general Yang-Mills theory with massive 
gauge bosons. Free massive spin 1 fields with color index a satisfy the Proca equation 

nA^,-d>^id,A:) + mlA^,=0. (88) 

Taking the divergence of H88|) shows that the field automatically satisfies the Lorentz 
condition df^Aa = 0. Again we may quantize the fields as four independent scalar fields 
which fulfill the wave equation 

{n + ml)A^,=0, (89) 

and the Lorentz condition shall be enforced on the Fock space like in the massless case by 
selecting a physical subspace which contains the three allowed polarizations. The definition 
of the gauge charge via does no longer lead to a nilpotent gauge charge due to the 
mass of the gauge field. This can be restored by introducing scalar bosonic fields with 
the same mass as the corresponding gauge field 

(□ + m2)$„ = 0. (90) 

These fields are called Higgs ghosts or would-be Goldstone bosons, and they are unphysical. 
We demand that the free fields then satisfy the commutation relations 

mx),AUy)]=iSab9^''D^^{x-y), (91) 

[Ua{x),Ub{y)] = -i^abDmaix - y), (92) 

[^a{x),My)] = -iSabD„,Ax-y), (93) 

where — y) are the usual Pauli-Jordan distributions for massive fields 

Dm{x) = I d'k S{e - m2)sgn(fc0)e-'*=^ (94) 

Then we can define a nilpotent gauge charge by 

Q := J d^x (9^A^(x) + ma^a{x)) do u{x) (95) 
Gauge variations are then obtained from the commutators 

[Q,A^]=td>'Ua, {Q,Ua} = 0, 
{Q, Ua} = -i{df,AI^ + ma<^a), [Q , ^a] = irUaUa. (96) 

The coupling H63|) which we obtained in the massless case 



Tiix) = igfabcU ■■ A^aix)A,bix)F^^'{x) : - : A^aix)ubix)d''ucix) :} (97) 
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is not gauge invariant anymore due to the additional massterms in the gauge variations 
(|96() . But gauge invariance can be restored at least at first order if we add to Ti a Higgs 
ghost couphng 



for the case where ah masses are nonvanishing. But this is not the full story, because gauge 
invariance can be shown to break down again at second order under certain circumstances. 
As a consequence, a complete analysis of the situation is necessary also at higher orders. 
In order to illustrate the implications of gauge invariance, we present here the result of 
such an analysis for the case where we have three different colors, and we assume that 
at least one of them is massive (because otherwise we would be dealing with a massless 
SU(2) Yang-Mills theory). 

It is found that there exist two minimal scenarios, which are gauge invariant also at 
second order, whereas the coupling fate is equal to the completely antisymmetric tensor 
Eabc for both cases: 

One possibility is given by the case where two gauge boson are equally massive, and 
the third one is massless. The correct coupling structure can be obtained from ()98() by 
performing a limit ms — > for mi = m2- The corresponding massless Higgs ghost does 
not appear anymore in the gauge charge Q, and becomes a massless physical field. The 
second possibility is given by the case where all three gauge bosons are massive, and all 
masses are then equal. But it turns also out that it is necessary to introduce additional 
fields in order to save gauge invariance. Indeed, it is already sufficient to add one physical 
scalar to the theory (i.e. this is the minimal solution): A Higgs boson. Gauge invariance 
fixes the structure of the coupling completely, but the mass of the Higgs remains a free 
parameter in the theory. The two scenarios correspond to the classical picture of the two 
possible types of spontaneous symmetry breaking for a SU(2) Yang-Mills theory: In the 
first case, one couples a real S0(3) field triplet to the gauge fields. Two degrees of freedom 
of this triplet are 'eaten up' by the gauge bosons, which become massive. In the second 
case, one couples a complex doublet (equivalent to four real fields) to the gauge bosons. 
Three degrees of freedom get absorbed by the gauge fields, which become massive, and 
one degree of freedom remains as a physical Higgs field. These two mechanisms are well 
described in many textbooks (see e.g. [231, chapter 8.3). A detailed discussion in the 
present approach of theories involving massive gauge fields like the Z- and T^^-bosons 
can be found in j251 1^6 . 27.. 28, 29,- The interesting point of the free field approach is 
given by the fact that it reverses the chain of arguments usually presented in the literature. 
Fields do not become massive due to the coupling to a Higgs field, but the existence of a 
Higgs field becomes necessary due to the mass of the gauge bosons and as a consequence 
of perturbative gauge invariance. 

8 Quantum gravity (as massless spin 2 gauge theory) 

For the following discussion, we should always keep in mind that quantum gravity is a 
nonrenormalizable theory, and not really understood to all orders in perturbation theory. 





(98) 
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But at least it is possible to discuss gauge invariance of the lowest orders PUII^ . 

The general theory of relativity can be derived from the Einstein-Hilbert Lagrangian 

J^EH = -^V^R (99) 

where R = g^^ R^u is the Ricci scalar, g is the determinant of the metric tensor g^^^ and 
K? = 327rG, where G is Newton's gravitational constant. It is convenient to introduce 
Goldberg variables ^| 

= V^g"". (100) 

Since we are working on a perturbative level, we expand the inverse metric g'^'^ and assume 
that we have an asymptotically flat geometry 

Here rj^^ = diag(+l, —1, —1, —1) is the metric of Minkowski spacetime. We used the more 
common symbol g^'^ for the flat space metric in the previous sections, but use rj^" in the 
following in order to point out the difference between ry and g more clearly. All tensor 
indices will be raised and lowered with •q^'^ . The quantity h^'^ is a symmetric second rank 
tensor field, which describes gravitons after quantization. For notational convenience, we 
do not care if the indices are up or down, since the summation convention always makes 
clear what is meant. Formally ()99() can be expanded as an infinite power series in n 

oo 

i=o 

The lowest order term £^1 is quadratic in h^'^{x) and defines the wave equation of the 
free asymptotic fields. The linearized Euler-Lagrange equations of motion for h^^{x) are 

Uh^'^ix) - ^r7^^n/i(x) - h^^P^^{x) - h'^/^'^ix) = (103) 

where h{x) = h^{x), and the commas denote partial derivatives with respect to the corre- 
sponding index. This equation is invariant under a classical gauge transformations of the 
form 

h^'^ ht"^ + u''^" + u"^'' - ri^'^uPp (104) 
where is a vector field which satisfies the wave equation 

Uu^{x) = 0. (105) 

As a consequence, quantum gauge invariance for gravity will be formulated by the help of 
(fermionic) vector ghosts. This gauge transformation corresponds to the general covariance 
in its linearized form of the metric tensor g^jy{x). The corresponding gauge condition, 
compatible with (|104|) is the Hilbert-gauge, which obviously plays a similar role as the 
Lorentz gauge condition for spin 1 fields 



(106) 
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Then the dynamical equation for the graviton field h^'^ reduces to the wave equation 

ahf^'ix) = 0. (107) 
The first order term with respect to k, C^^]j gives the trilinear self-couphng of the gravitons 

c'^eI = \ h^" {hfh'f - i h^h^ + 2 hjh^: + h^h^S - 2 hjhj) . (108) 

There exists many alternative derivations of this result ()108() . starting from massless 
tensor fields and requiring consistency with gauge invariance in some sense. A similar 
point of view is the one of Ogievetsky and Polubarinov They analyzed spin 2 theories 
by working with a generalized Hilbert-gauge condition to exclude the spin one part from 
the outset. They imposed an invariance under infinitesimal gauge transformations of the 
form 

dh^"" = d^'u" + d'^u'' + r]^"'^au'' (109) 

and obtained Einstein's theory at the end. Instead Wyss [HII considered the coupling to 
matter. Then the self-coupling of the tensor- field (|108() is necessary for consistency. Wald 
|35j derived a divergence identity which is equivalent to an infinitesimal gauge invariance 
of the theory. Einstein's theory is the only non-trivial solution of this identity. In quantum 
theory the problem was studied by Boulware and Deser [301 • These authors require Ward 
identities associated with the graviton propagator to implement gauge invariance. All 
authors get Einstein's theory as the unique classical limit if the theory is purely spin two 
(without a spin one admixture). 

The free asymptotic field h'^'^ is a symmetric tensor field of rank two and and are 
ghost and anti-ghost fields on the background of Minkowski spacetime in the following. A 
symmetric tensor field has ten degrees of freedom, which are more than the five indepen- 
dent components of a massive spin 2 field or the two degrees of freedom of the massless 
graviton field. The additional degrees of freedom can be eliminated on the classical level 
by imposing further conditions, but we will now quantize the graviton field and introduce 
the gauge charge operator Q. The physical states in the full Fock space which contains 
eight unphysical polarizations of the graviton and ghosts can then be characterized di- 
rectly by the help of the gauge charge as in the case of a spin 1 field (see Appendix C, and 

jsniHiEni). 

The tensor field h^'^{x) can be quantized as a massless field as follows 

[h^^{x),h>'-{y)] = -ib-f^f^-'Doix-y) (110) 

where Dq{x — y) is again the massless Pauli- Jordan distribution and the tensor b°'^'^'^ is 
constructed from the Minkowski metric rj^'^ via 

b'"^^"' = i (r/"^/^ + r/"'^/^ - 7?°^7?^'^) . (Ill) 
An explicit representation of the field h^^^{x) on the Fock space is given by 

h^'^ix) = (27r)-3/2 ( ^ (a°^(fc)e-''=" + a'^^ [k)^ e+'^A . (112) 
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Here we have u; = |A:| as usual and a°'^{k), a^^^k)^ are annihilation and creation operators 
on a bosonic Fock-space. One finds that they have the following commutation relations 

[a'"^{k),a^"'{k')^ =b"^^''S^'^\k-k'), (113) 

and all other commutators vanish. 

In analogy to the spin 1 theories like QCD we may introduce a gauge charge operator 

by 

Q:= y d^xh'^'^{x)^f3doUa{x). (114) 

For the construction of the physical subspace and in order to prove the unitarity of the 
S'-matrix we must have a nilpotent operator Q. Therefore we have to quantize the ghost 
fields with anticommutators according to 

{u^{x),u''{y)}=iv^''Do{x-y). (115) 

All asymptotic fields fulfill the wave equation 

□ /i/^'^(x) = □u"(x) = □u^(x) = 0. (116) 

We obtain the following commutators of the fundamental fields: 

[Q, h] = iu^^, (118) 

{Q,n^} = ^/l(t^ (119) 

{Q,u^} = 0. (120) 

From piTj) we immediately get 

[Q,h':,n = o. (121) 

The result ()117() agrees with the infinitesimal gauge transformations of the Goldberg vari- 
ables, so that our quantization and choice of Q corresponds to the classical framework. 

Again, first order gauge invariance means that [Q, Ti] is a divergence in the sense of 
vector analysis, i.e. 

[Q,Ti{x)]=id^T^^^{x), (122) 
and definition of the n-th order gauge invariance then reads 

n Q 

[Q,Tn] = [Q,Tn\ =«^— ^r^/,(xi,...,Xi,...,x„). (123) 
1=1 

Like in the case of QCD, we are forced to add a ghost part to the self-coupling term of 
the gravitons. The first order graviton coupling according to ()1U8() is 

= ^ h"^ (hfhf - \ h^h^ + 2 hT^^h^: + h^h':^ - 2 hT^^h'^^) . (124) 
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By adding a physically irrelevant divergence to T{*, we obtain a more compact expression 
for the graviton interaction 



The ghost coupling tm'ns out to be the one first suggested by Kugo and Ojima [301, namely 



It is a nice detail that the four graviton vertex which follows from the second order term in 
()102() is also proliferated by gauge invariance at second order. It is therefore quite probable 
that also all higher vertex couplings which appear in ()102j) are proliferated by quantum 
gauge invariance. It is also possible to derive the couplings H124() . H126|) from perturbative 
gauge invariance similarly as it was done in sect. 6 for Yang-Mills theories. For this and 
further reading we refer to the recent monograph j22j . 

We finally point out that it would be interesting to analyze the Higgs mechanism for 
massive gravity [IJI in the framework presented above, also in connection with the dark 
matter problem. A first step in this direction is presented in |42j . 

9 Appendices 

Appendix A: Commutator of the gauge charge Q with gauge fields 

First we derive the distributional identity (|16|1. The Pauli- Jordan distribution has an 
integral representation 




(125) 




(126) 




(127) 



Using the identity 




(128) 



we obtain (mind the sgn(/c'^)) 



doD{x) 



i f d^k 



[Sik' - \k\) 



5(A:° + |fc|)(-ife' 



(27r)3 J 2\k0\ 



2(27r)3 y V" 
Restricting this result to = 0, we get the desired result ()16j) 




(129) 




(130) 



In a completely analogous way, one derives 



diD{x%o=o = 0, VZ)(x)Uo=o = 0. 



(131) 



9 APPENDICES 



21 



Note that we always consider the weU-defined differentiated distribution first, which gets 
then restricted to a subset of its support. 

As an example, we investigate now the commutator [Q, A^{y)], and we omit the trivial 
color index. The commutator is given explicitly by 

[Q,A^{y)] = [ [ d^xd,A''{x)dSu{x),A^{y)]=i [ d'x d^^D{x - y) u{x). (132) 



Here, we made use of the freedom to choose any constant value for x^, and therefore we 
set x^ = y^, such that we have x^ — y^ = and we can apply ((TBI I131|) in the sequel. 
For /i = 0, we have 

'3^ ax \ 



[Q,A^{y)]=i / d-'xd^^D{x-y)d^^u{x) 



d^x5^^\x-y)d^u{x) = idou{y), (133) 

where we have used that the double timelike derivative of D vanishes on the integration 
domain according to (|131jl . The result for the commutator of Q with the spacelike com- 
ponents of A is also obtained by using (|161 1131|l and by shifting the gradient acting of the 
Pauli- Jordan distribution by partial integration on the ghost field. 

Appendix B: The Becchi-Rouet-Stora-Tyutin transformation and its free 
field version 

The gluon vector potential can be represented by the traceless hermitian 3x3 standard 
Gell-Mann matrices A", a = 1, ...8 

^m = EA".Y=:^^Y- (134) 

a=l 

The A's satisfy the commutation and normalization relations 



A'' A'' 
Y' Y 



i/afecy, tr (A'^A") = 2,5,6, (135) 



and the numerical values of the structure constants fabc = —fbac = — fact can be found in 
numerous QCD textbooks. Since we are working with a fixed matrix representation, we 
do not care whether the color indices are upper or lower indices. 

The natural generalization of the QED Lagrangian to the Lagrangian of purely gluonic 
QCD is 

1 u a u 

C-giuon = — -tr G^jyG'^'^ = — -G^jyG'^'^, (136) 

with 

G^. = d^A, - d,A^ - ig[A^, A,] (137) 
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or, using the first relation of (|135() 

G% = d^Al - + gfahcA^Al. (138) 

It is an important detail that we are working with interacting classical fields here, therefore 
the field strength tensor contains a term proportional to the coupling constant in contrast 
to the free field tensor F^il^^ = d^Al^^^ — dyAfl'^^ used in this paper. Cgiuon is invariant 
under classical local gauge transformations 

A^{x) ^ U{x)A,,{x)U-\x) + -U{x)d^U-\x), (139) 

g 

where U{x) e 5C/(3). 

We extract now the first order gluon coupling from the Lagrangian. The Lagrangian 

Cgiuon = -\[d^.Al - d^Al + gfabcAlAlWAl - d''A^,+gfab'c'A'^,A''^,], (140) 

contains obviously the free field part (this terminology is not really correct, since we are 
dealing with interacting fields here) 

= -li^.Al - d^A-WA: - a^^] (141) 
and the first order interaction part is given by 

Clfuon = -\[d,Al - d,Al][gf,,,,A>^,,A:,] - \[g fatcA^AlW A^ - d'^A^J 
= -^f^,^AlAl[d^A: - d^A^J = -3-f,,^AlAl[d^A- - A^^] 

= gfabcA'lAld-A'::. (142) 

From this term follows the first order interaction ansatz = i : >C^"*q„ : H53() . 

Since we are working in Feynman gauge, we add the corresponding gauge fixing term 
Cgf to the Lagrangian. Additionally, we add a ghost term which leads to the ghost 
interaction (|59j) . The total Lagrangian then reads 

C-QCD = Cgiuon + Cgf + Cghost 



= Cgiuon - ^{d^Ai^f + df^uid^Ua " S/afec^feV)- (1^3) 

The classical ghosts are anticommuting Grassmann numbers, i.e. u'^ = v? = 0, uu = —uu. 
The BRST transformation is defined by 

dA"^ = i\{d^ua - gfabcUbA^c), (144) 
6ua = -i\df,Ai^, (145) 



9 APPENDICES 



23 



5Ua = ^XfabcUbUc, (146) 

where A is a space-time independent anticommuting Grassmann variable. The special 
property of the BRST transformation is the fact that the actions 

Sgluon — j d'^X ^gluoni Sgf + Sghost — J d'^^ i^gf ~l~ ^ghost) (147) 

and Stotai = Sgluon + Sgf + Sghost are all invariant under the transformation: 

SSgluon = 0, 6iSgf + Sghost) = 0. (l48) 

The similarity of free quantum gauge transformation introduced in this paper to the 
BRST transformation is obvious. One important difference is the absence of interaction 
terms ~ g. Furthermore, the free quantum gauge transformation is a transformation of 
free quantum fields, whereas the BRST transformation is a transformation of classical 
fields, which enter in path integrals when the theory is quantized. Finally, the free gauge 
transformation leaves the T„'s invariant up to divergences, whereas the BRST transforma- 
tion is a symmetry of the full QCD Lagrangian. How the two symmetries are intertwined 
perturbatively is explained in |^. A more rigorous axiomatic approach is discussed in 

113 uni. 



Appendix C: Some technical remarks concerning the gauge charge Q 

Using the Leibnitz rule for graded algebras gives for the gauge charge for massless spin 1 
fields 

= liQ.Q} = \ j d\dMx){d., Ua{x),Q} 

xo=const. 

d'x[d,A'',{x),Q]f,,Ua{x) = , (149) 

xo=const. 

i.e. Q is nilpotent. This basic property of Q, which holds also in the massive case, and the 
so-called Krein structure on the Fock-Hilbert space |441 145j allows to prove unitarity of 
the 5-matrix on the physical Hilbert space J-phys, which is a subspace of the Fock-Hilbert 
space J- containing also the unphysical ghosts and unphysical degrees of freedom of the 
vector field 

The physical Fock space can be expressed by the kernel and the range of Q [46| I44j . 
namely 

Fphys = kerQ/ranQ = ker{Q,Q'^] . (150) 

Again, this characterization of the physical space holds also for the massive spin 1 case, and 
for the gauge charge (|114)) used for gravity. The Krein structure is defined by introducing 
a conjugation K 

ao{k)^ = -ao{k)\ aj{k)^ = aj{k)^ , (151) 
so that A^ = A^^ and on the ghost sector 




C2{kf = Ci{k)\ Cl(fc)^ = C2(^)^ 



(152) 
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so that u = u is i^-selfadjoint and u = —u. Then Q is densely defined on the Fock- 
Hilbert space and becomes i^-symmetric Q C . Roughly speaking, the K-conjugation 
is the natural generalization of the usual hermitian conjugation to the full (unphysical) 
Fock space. 

A calculation shows that the anticommutator in is essentially the number oper- 

ator for unphysical particles 

{Q\Q} = 2 J d3A:P[6|(fe)6i(fe) + 6^(fc)62(fc) + cl(fc)ci(fc) + 4(^)c2(^)], (153) 

with 

6i,2 = (ay ± ao)/\/2 , aii=kjaj/\k\, (154) 

which implies 

The nilpotency of Q allows for standard homological notions 0^1 • Consider the field 
algebra J- consisting of the polynomials in the (smeared) gauge and ghost fields and their 
Wick powers. Defining a gauge variation for a Wick monomial F according to 

dQF=^QF-{-ir^FQ, (155) 

where np is the number of ghost fields in F, Q becomes a differential operator in the sense 
of homological algebra, and we have 

4 = 0^ {Q, [Q, Ft]} = [Q, {Q, Ff}] = 0, (156) 

where Fb is a bosonic and Ff a fermionic operator and dQ{FG) = {dQF)G FdqG. 
For example, we get 

dg : Af^aUbd^Uc ■■=: [Q, Af,a]ubd^Uc : + : A^alQ, UbWuc : - : Af,aUb{Q, d^'uc} : . (157) 
If F = dqG, then F is called a coboundary. The term 

fabc ■ duA'^UbUc ■■= -dgifabc ■ UaUbUc :) . (158) 

is such a coboundary, which does not contribute to the physical QCD S-matrix. 



Appendix D: Gauge invar iance at n-th order 

A short discussion is given here which explains how the simple condition of gauge invariance 
at n-th order in the coupling constant (|46j) emerges as a generalization from the first order 
condition (|lHl . 

A thorough mathematical treatment of the subject would involve a discussion of op- 
erator valued distributions, which is not given here for the sake of brevity. Nevertheless, 
all mathematical steps presented in the following can be put on a sound mathematical 
basis. It is useful to remind some basic facts: Free quantum fields like the free scalar field 
f{x) @ are operator valued distributions on a Fock-Hilbert space, i.e. an operator '^{g) 
is obtained after smearing out the field operator (p{x) with a test function g{x). This can 
be written formally as 

ip{g) = [ d\ip{x)g{x), (159) 
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where g G 5(R'^) is in the Schwartz space of infinitely differentiable and rapidly decreasing 
test functions. Furthermore, the tensor product ip{x)ip{y) of free fields (and normally 
ordered products of free fields) is also an operator valued distribution, whereas the local 
product ip{x)'^{x) is only defined after normal ordering. The 'value' of a distribution in a 
single point makes no sense in general. It is sometimes convenient to treat distributions 
like ordinary functions, as we will do it below, because in many cases the simplified insight 
gained from such a simplification allows to construct a full mathematical proof at a later 
stage. 

As mentioned in the paper, the r„(a;i, ...x„) are well-defined time-ordered products of 
the first order coupling Ti(x), and they are expressed in terms of Wick monomials of free 
fields. The construction of the T„ requires some care: If the arguments {xi, ■■■Xn) are all 
time-ordered, i.e. if we have 

xl>xl> ... > 4, (160) 

then Tn is simply given by 

Tn{xi,...Xn) =Ti{xi)Ti{x2)...Ti{Xn) (161) 

According to the definition the T„(xi, ...Xn) can be considered symmetric in xi, ...Xn. 
Using this fact allows us in principle to obtain the operator-valued distribution r„ induc- 
tively everywhere except for the 'complete diagonal' A„ = {xi = ... = Xn} 03- The 
construction is inductive since in the construction of T„ subdiagrams with lower order 
appear. If T„ were a C-number distribution, we could make it a well-defined distribution 
for all xi, ...Xn by extending the distribution from 7^^"/A„ to TZ'^^. In the case of free field 
operators, the problem can be reduced to a C-number problem by the Wick expansion of 
the operator-valued distributions. The extension T{xi, ...Xn) is, of course, not unique: it 
is ambiguous up to distributions with local support A„. This ambiguity can be further 
reduced by the help of symmetries (in particular gauge invariance) and power counting 
theory, and it is this local ambiguity which shows up as ultraviolet divergences in Feynman 
diagram calculations. 

The concrete inductive construction of the T„ is discussed in detail in a famous paper 
of Epstein and Glaser |3H| for scalar field theories, and explained in a more pedagogical 
way in [IHI- For many practical calculations, it is advantageous to work in momentum 
space since the distributions, i.e. the Green's functions which occur in the calculations, 
behave much smoother in p-space than in x-space. The work of Epstein and Glaser is 
based on a treatment in real space, whereas in practical calculations are performed 
in momentum space. An advantage of the real space formulation is also that it allows 
to formulate a consistent renormalization theory on a curved physical background j5,Sj . 
Interesting topics like scale invariance, renormalization and the renormalization group can 
also be treated in real space jS21IS3I^- 

As a specific example, we consider the second order contribution to the S-matrix which 
is given by 

S2 = ^J d'xd^yT[T,{x)T,{y)] = j d^x d''yT2{x,y). (162) 
We consider first the simple product Ti{x)Ti{y) without time ordering. Then we have 



[Q,T,{x)T,{y)] = [Q,ri(x)]ri(y) +ri(x)[Q,ri(y)] 
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= i5^r{)^(x)Ti(y) + iTi(x)9^r()^(y). (163) 

Considering the time-ordered product T[Ti{x)Ti(y)], we have to distinguish three cases. 
In the first case, we have > y^, and therefore T[Ti{x)Ti{y)] = Ti{x)Ti{y). The second, 
a bit less trivial case is given for x^ < y^ ■ Then we have T\Ti{x)Ti{yy\ = Ti{y)Ti{x), and 

[Q,T{Ti{x)Tiiy))] = [Q,Ti{y)]Ti{x) + Ti{y)[Q,Ti{x)] 

= id^T^Jy)Tiix) + iTi{y)d^T^,^{x) 



= ir[9^r{)^(x)ri(y)] + iT[T^{x)d,T';^^{y)]. (164) 

Finally, we may have x^ = y^. If {x — y) ^ 0, i.e. if x and y are spacelike separated, then 
we may perform 'a little' Lorentz transformation such that x^ > y^, and the discussion 
above applies. If the special case that x = y, gauge invariance is not trivially fulfilled. This 
is not astonishing, since this is the case on the diagonal A2 = {x = y} where perturbation 
theory may fail to calculate time-ordered products of the first order interaction. Requiring 
that gauge invariance holds also for x = y leads to Ward-Takahashi identities in QED, 
and to the so-called Slavnov- Taylor identities for QCD. 

The simplest example for such an identity is obtained from the vacuum polarization 
contribution to the S marix at second order, which can be written in the form 

T^^{x,y)=:A^{x)t^^''{x-y)A,iy) :, (165) 

where t^'^{x — y) = t'^^{x — y) is a C-number distribution. Operator gauge invariance 
requires 

[Q,rf (x,y)] = i9^n(x)t^^(x - y)A,{y) + iA^(x)t'^^(x - y)9.n(y) 



= id-^[n{x)t'^''{x - y)A,{y)] + i9^[^^(x)t^^(x - y)u{y)], (166) 

and therefore 

a^t^^(x-y) = a,^t^'^(x-y) = 0. (167) 
In momentum space, we get the C-number identity 

k^i^'^ik) = kJ^"'{k) = 0, (168) 

hence the vacuum polarization tensor is transversal. 

Appendix E: The adiabatic limit 

The perturbative expression is problematic, because the time-ordered products r„ 
are operator valued distributions after regularization, and they have to be smeared out by 
test functions. In order to be more precise in a mathematical sense, we introduce a test 
function go{x) £ 5(R^) with go{0) = 1 and replace expression by 

5 = 1 + ^— / d^xi . . .d'^XnTnixi, . . .Xn)goixi)...goixn)- (169) 
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Here, go acts as an infrared regulator, which switches off the long range part of the 
interaction in theories where massless fields are involved. E.g., in QED the emission 
of soft photons is switched by go, and as long as we do not perform a so-called adiabatic 
limit go ^ 1, the matrix elements of the S-matrix remain finite. One possibility to perform 
the adiabatic limit is by scaling the switching function go{x), i.e. one replaces go{x) by 
g{x) = goi^x) and performs the limit e — > 0, such that g approaches everywhere the value 
1. If the S-matrix is modified by a gauge transformation, operators which are divergences 
are added to the n-th order term r„. Such a contribution can be written as 

J (fxi...(fxnd''^^0-^-{xi, ...,xi, ...Xn)g{xi)...g{xi)...g{xn) 

= - j (fxi...d'^XnO-^'-{xi, ...,xi, ...Xn)g{xi)...dl'g{xi)...g{xn). (170) 

In the adiabatic limit, the gradient d^'-g{xi) vanishes. Unfortunately, this property of the 
scaling limit does not guarantee that the whole term (|170|) vanishes (see also [SUES])- 

The infrared problem is not really understood in QCD, and all proofs of unitarity which 
exist in the literature have to be taken with a grain of salt, because they are avoiding the 
discussion of infrared problems somehow. 

A thorough perturbative approach to the construction of the local algebras of observ- 
ables which avoids the adiabatic limit is given for QED in and (under the assumption 
that there are no anomalies) for non-Abelian gauge theories in |19j . 
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